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ABSTRACT 


This  paper  discusses  recent  work  on  the  development  of  analysis 
procedures  for  obtaining  drift  and  deformation  measured  from  sequential 
visual  imagery  of  sea  ice  that  is  located  far  from  land.  In  particular 
for  LANDSAT  images  far  from  land  a semi  automatic  procedure  for  trans- 
ferring the  location  coordinates  of  a common  set  of  ice  features  from 
the  Earth  coordinate  system  of  one  image  to  another  is  discussed. 
Necessary  inputs  for  the  transfer  are  the  location  coordinates  (lati- 
tude and  longitude)  of  the  center  of  each  image  and  the  location  of 
two  arbitrary  points  on  a known  line  of  longitude;  all  this  information 
is  available  from  LANDSAT,  although  with  some  error.  Errors  in  the 
transferral  technique  are  examined  using  imagery  over  land  and  are 
found  to  be  dominated  by  deviations  (as  large  as  8 km)  in  the  actual 
position  of  the  center  of  the  image  from  its  stated  position.  The 
errors  on  the  average  are,  however,  less  than  typical  one  day  ice  drift 
distances.  The  LANDSAT  image  location  errors  also  introduce  uncertain- 
ties in  the  orientation  of  the  coordinate  systems  after  transfer. 

These  errors  will  produce  spurious  apparent  strains  if  velocities  are 
estimated  by  simply  taking  position  differences.  This  subtle  effect 
may  also  occur  if  the  coordinate  system  of  a given  floe  is  used  as  the 
common  coordinate  system.  A least-squares  strain  program  utilizing 
polar  coordinates,  which  eliminates  such  spurious  effects,  is  dis- 
cussed. 

With  regard  to  measuring  strain  from  sea  ice  aerial  imagery  with- 
out ground  control,  errors  in  such  measurements  are  examined  using 
uncorrected  photographs  obtained  during  April  1975  by  MSA  using  a 
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RC-8  camera  over  the  main  AIDJEX  camp.  During  the  time  of  the  photo- 
graphic flights  a '^7  kilometer  "calibration"  line  vas  monitored  on  the 
ground.  The  errors  in  using  such  uncorrected  imagery  and  using  common 
undeformed  ice  floes  to  establish  a common  scale  are  found  to  be  of  the 
order  of  1%  vhereas  typiceJ.  maximum  differential  motions  are  as  large 
as  5^. 

INTRODUCTION 

One  use  of  the  satellite  and  aerial  imagery  of  considerable 
importance  for  modelling  studies  is  the  estimation  of  sea  ice  drift  £uid 
deformation  by  using  sequentieil  overlapping  Images.  Such  a use  of 
LANDSAT  imagery  near  land  has  been  discussed  by  a number  of  authors  (eg. 
Crovder  et.  al. , (197^);  Shapiro  and  Burns,  1975).  In  these  papers 
the  basic  procedure  has  been  to  identify  common  features  on  sequential 
images  so  that  relative  distance  changes  can  then  be  measured.  If 
some  points  on  the  image  of  interest  are  on  land  then  drift  rates 
in  addition  to  deformation  rates  may  be  estimated  since  the  land  points 
do  not  move. 

However,  when  LANDSAT  Images  far  from  shore  are  used  for  analysis, 
the  problem  of  obtaining  both  drift  and  strain  rates  becomes  consider- 
ably more  difficult.  The  principal  difficulty  with  drift  rate  estima- 
tion is  that  no  immovable  land  points  are  available  for  "calibration" 
and  consequently  one  has  to  use  the  coordinates  of  the  centers  of  the 
images  as  well  as  available  longitude  marks  to  fix  the  image  locations. 
In  principal  the  problem  of  drift  then  becomes  one  of  first  projecting 
common  ice  featiires  (on  sequential  images)  onto  the  spheroid  and  then 
back  onto  a common  tangent  plane.  A somewhat  more  useful  procedure, 
which  we  will  discuss  in  this  paper,  is  to  transfer  the  coordinates  of 
both  images  (considered  to  be  tangent  planes)  onto  a common  inter- 
mediately located,  tangent  plane.  In  practice,  such  a procedure  simply 
consists  of  appropriately  translating  and  rotating  the  polar  coordinates 
of  all  points  in  sequential  images  to  the  coordinate  system  of  the 
first  image. 

The  estimation  of  the  strain  rate,  on  the  other  hand,  is  depen- 
dent in  a more  subtle  way  on  the  transfer  of  coordinates  from  one 
Imaje  to  another.  The  key  problem  here  is  caused  by  spurious  rotations 
of  the  coordinate  system  that  cm  be  induced  by  errors  in  longitude 
lines.  Such  spurious  rotations  can  also  be  induced  by  using  a common 
rotating  floe  to  establish  a coomon  coordinate  system.  The  basic 
problem  is  that  if  velocities  are  estimated  by  merely  subtracting 
location  coordinates  in  a rectangular*  coordinate  system  then  such 
spurious  rotations  will  Induce  concomitant  errors  in  the  diagonal 
compunents  of  the  strain  rate  tensor.  Such  undesirable  effects  cam 
be  avoided  by  using  an  appropriate  least-squares  strain  program  in 
polar  coordinates,  an  example  of  which  will  be  discussed  later. 

In  the  following  we  will  discuss  procedures  used  to  estimate  the 
drift  and  deformation  of  sea  ice  from  lANDSAT  images  of  areas  located 
far  from  land.  We  then  briefly  examine  some  results  obtained  from 
uncorrected  aerial  imagery  photographed  from  91^^  m (30,000  ft)  altitude. 
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DRIFT  AND  DKFORMATION  FROM  LANDSAT  IMAGERY 


Drift  Rate  Estimation 

In  principle,  if  the  center  locations  of  two  sequential  over- 
lapping LANDSAT  images  are  known,  then  drift  rates  ma7  he  estimated 
by  identification  of  ice  features  common  to  both  images.  One  way  to 
do  this  is  to  transfer  the  coordinates  from  the  image  planes  onto  the 
spheroid,  and  then  onto  a common  plane.  Such  a procedure  is,  however, 
relatively  complex.  A simpler  procedure,  adequate  for  LANDSAT  imagery 
is  to  assume  that  both  images  are  effectively  co-planar  so  that  the 
transfer  of  coordinates  becomes  a simple  translation  and  rotation  in  a 
given  plane. 

In  particular  it  is  convenient  to  asstuae  that  both  images  are  in 
the  plane  whose  perpendicular  is  given  by  the  cross  product  of  the  two 
north  vectors  tangent  to  the  earth  at  latitude  6 , and  parallel  to  the 
longitude  lines  ^ and  The  angle  between  these  two  vectors  defines 

the  angle  between  the  x-axes  of  the  two  image  coordinate  systems  for 
image  centers  on  the  same  latitude  line.  For  simplicity  we  also  take 
this  angle  to  be  the  same  if  the  second  image  is  on  both  a different 
longitude  line  and  a different  latitude  than  the  first , an  assumption 
inducing  errors  of  order  d^/R^.  By  taking  an  inner  product,  the  angle 
between  these  two  vectors  y,  is  given  by 

—1  (2  2 

y = cos  |cos  6 + sin  6 cos 
For  small  values  of  this  expression  reduces  to 

li^cos”^  |cos  sin  ((1]^=  sin  ip 


With  the  above  assumptions  the  geometry  for  the  coordinate  trans- 
fer is  given  in  Figure  1,  where  points  1 and  2 are  the  centers  of 
Images  1 and  2 respectively.  The  angle  of  rotation  between  the 
coordinate  systems  is  given  by  y and  the  translation  distance  is  defined 
by  the  chord  distances  Dj^  and  Dp  along  the  latitude  and  longitude  lines. 
The  values  of  Dj^  and  Dp  can  be  calculated  from  the  expressions: 


2 sin 


LAT 


°2  = 


2 sin 

[e*-el 


'e*-e~ 

2 . 


LONG 
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where  Dlahi  is  the  dlstemce  on  the  International  Spheroid  between  longi- 
tude <f'  and  measured  along  latitude  line  0,  and  i>XiQNG  inter- 

national Spheroid  distance  between  latitudes  0 ' and  0 measured  along 
longitude  line  (|i*.  Approximate  expressions  for  these  distances  are 
given  on  page  ll87  of  Bowditch  (1958). 

In  order  to  give  an  idea  of  the  errors  induced  by  the  approxima- 
tions in  the  coordinate  transfer  procedure,  we  have  compared  the  dis- 
tance and  angles  between  a number  of  points  obtained  using  our  LANDSAT 
projection  procedure  for  both  a sphere  (referred  to  as  the  LANDSAT 
SPHERE)  and  the  International  Spheroid  (referred  to  as  the  LANDSAT 
spheroid)  with  the  results  obtained  by  using  a planar  projection  onto 
a fixed  pleuie  tangent  to  point  1.  For  the  spherical  cases,  we  con- 
sidered point  1 to  be  at  lat  80°,  long  0°  and  have  also  considered  9 
other  points,  3 located  along  lat  80°  at  increasing  distances  towaj-d 
the  west,  3 located  along  long  0°  at  increasing  distances  toward  the 
north  and  3 located  at  increasing  distances  toward  the  northwest.  The 
results  are  summarized  in  Table  I with  positive  angles  being  measured 
in  a clockwise  direction  from  the  x (north)  axis.  Point  1 is  taken  to 
be  the  origin.  In  the  spherical  calculations  a radius  equal  to  the 
average  radius  of  the  International  Spheroid  was  used.  The  great 
circle  distance  was  computed  using  the  spherical  surface.  In  general. 
Table  I shows  that  there  is  very  little  difference  between  the  various 
projections  with  the  maximum  difference  being  about  1 km  at  600  km. 

In  actual  practice  the  maximum  transfer  distances  are  rscrely  greater 
than  250  km  so  that  differences  in  the  various  projections  are  hardly 
significant. 

In  order  to  test  the  transferral  technique  in  actueil  practice  we 
utilized  a number  of  overlapping  LANDSAT  images  of  Alaska.  From  pairs 
of  such  imEiges  ccmimon  land  features  were  identified  and  their  coordi- 
nates were  transferred  from  the  second  image  to  the  coordinate  system 
of  the  first  image.  Because  the  points  were  "immovable”,  apparent 
motions  could  be  identified  as  errors. 

After  examination  of  a ntjmber  of  such  image  pairs  it  became 
apparent  that  the  dominant  source  of  error  in  the  coordinate  transfer 
was  caused  by  deviations  between  the  actual  location  of  the  center  of 
the  image  and  its  stated  location.  On  the  averaige  these  positional 
Inadequacies  induced  errors  of  the  order  of  one  kilometer. 

The  actual  errors  obtained  from  a series  of  6 photo  pairs  are 
given  in  Table  II.  In  this  table  ms  errors  before  and  after  correct- 
ing the  location  of  the  center  position  are  given.  The  center  posi- 
tions were  corrected  by  overlaying  the  images  on  accxjrate  survey  maps 
with  a zoom  transfer  scope.  As  can  be  seen  from  Table  II,  using  the 
actual  center  position  Improves  the  translation  errors  by  about  half  a 
kilometer  on  the  average.  The  actual  errors,  however,  between  the  true 
and  stated  center  location  can  be  as  large  as  8 km.  This  is  illus- 
trated in  Table  III,  which  summarizes  the  center  position  errors  of 
17  photographs.  The  maximum  deviation  obsejrved  was  8.12  km  (in  the 
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Figiire  1.  GEOMETRY  FOR  "LANDSAT"  PLANAR  TRANSFERRAL 


TABLE  I.  COMPARISON  OF  PLANAR  AND 
LAJiDSAT  PROJECTIONS. 


Location  of  initial  pt.  80”N,  0°W. 
Radius  of  Reference  Sphere  was  6'^T1.299  km 


Great 

LANDSAT  SPHEROID 

Location  of 

Circle 

LANDSAT  SPHERE 

Planar  Projection 

Projection  using 

Final  Pt 

Distance 

Projection 

using  Sphere 

International 

Spheroid* 

Long 

Lat 

D1 stance 

Angle 

Distance 

Angle 

Distance 

Angle 

Distance 

0®W 

81.2“S 

I33.U39  lun 

0° 

133.1*36  loti 

0“ 

133.1*29  Km 

0° 

133.572  Km 

o^w 

82.6°N 

289. IIT 

0° 

289. 092 

0“ 

289.018 

0° 

289.1*11* 

85  “N 

558.998 

0° 

555.819 

0° 

555.290 

0° 

556.1*15 

6°W 

80  °N 

115.806 

-87.01*6 

115.801* 

-87.0U5 

115.799 

-87.01*1* 

115.932 

111'’!! 

80  °N 

269.681 

-83.106 

269.661 

-83.105 

269.600 

-83.105 

269.959 

25°W 

80  °N 

1*79.030 

-77.690 

1*78.917 

-77.681* 

^78.579 

-77.688 

1*79.1*1*6 

6°w 

8i.2°N 

172.110 

-36.308 

172.113 

-36.302 

172.089 

-36.310 

172.291* 

ll*°W 

82.6"N 

370.871* 

-32.1*13 

370.915 

-32.390 

370.693 

-32.1*11 

371.323 

25°w 

85  °N 

651.1*77 

-21.262 

651.759 

-21.152 

650.31*2 

-21.261 

652.1*59 

•The  parametric  latitude  values  of  the  points  used  for  the  spherical 
calculations  were  converted  to  geographical  latitude  values  for  use  with, 
the  Spheroid  program. 


TABLE  II.  RKS  TRAMSLATION 

ERRORS. 

Photo  Pair 

Before  Center 
AX(North) 

Position  Correction 
AY (East) 

After  Center 
AX( North) 

Position  Correction 
AY ( East ) 

1 

1.66  Km 

l.lU  km 

.37  km 

.11  Km 

2 

2.ll* 

M 

1.22 

.1*9 

3 

1.63 

M 

.72 

.15 

U 

2.1*3 

.53 

.13 

5 

.55 

.11 

.55 

.58 

6 

.17 

.25 

.20 

.21* 

Average 

1.1*3  + .89 

.1*5  + .37 

.60  + .35 

.28  + .20 
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TABLE  III.  SUMMARY  OF  CENTER  POSITION  ERRORS  FOR 
17  IMAGES  OVER  ALASKA. 


X( North) 

Y(East ) 

Max.  Deviation 

1.73  km 

8.12  km 

Min.  Deviation 

0.20 

0.19 

RMS  Deviation 

0.84 

3.50 

east-west  direction).  These  numbers  are  in  general  agreement  with 
the  study  by  Colvocoresses  (197^)  who  has  found  1 to  8 km  errors  in 
the  latitude  and  longitude  center  indicators  on  LANDSAT  photos. 

The  reason  that  such  center  errors  do  not  induce  greater  errors 
in  the  transfer  appears  to  be  due  to  the  consistency  of  the  offset 
errors;  that  is,  pairs  of  photos  will  have  similar  errors  in  both  the 
magnitude  and  the  direction  of  their  center  positions.  However,  there 
is  no  assurance  that  the  updating  of  orbital  parameters  will  not  be 
made  in  the  middle  of  a sequence. 

Deformation  Rate  Estimation 

A convenient  way  of  estimating  least-squares  strain  rates  and 
vorticity  has  been  to  note  that  the  x and  y ice  velocities  at  position 
i,  (u  , u^  ) are  related  to  the  strain  rate  and  vorticity  by  (ass\jm- 
*i  ^i 

ing  a linear  velocity  field) 


A complete  discussion  of  the  application  of  least-squares  equations 
of  this  type  to  sea  ice  deformation  is  given  in  Hibler  et  al.  (197^). 
This  approach  has  the  advantage  that  the  velocity  points  further 
apart  automatically  affect  the  strain-rate  estimates  more  than  do 
velocities  of  points  close  together.  Also  the  vorticity  is  computed 
in  a normal  least-squares  way. 
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A drawback  to  this  approach  is  that  if  velocities  are  estimated 
by  merely  subtracting  location  coordinates  in  a rectangular  coordi- 
nate system  errors  will  be  induced  in  the  strain  rate  by  spurious  ro- 
tations of  the  coordinate  system.  The  basic  problem  can  be  easily 
visualized  by  examining  Figure  2.  In  it  we  consider  three  points, 
which  are  stationary  relative  to  one  another  and  consequently  undergo 
zero  strain.  Using  the  three  relative  distances  between  these  points 
we  would  in  fact  find  a strain  rate  of  zero  for  all  components  of  the 
strain  rate  rensor.  However,  using  rectangular  coordinates  and 
equations  (l)  we  would  obtain,  for  a clockwise  rotation  of  6 radians 
per  unit  time,  the  strain  rates  and  vorticity: 


"12  = ° 


^11  = "22  = 1 


w = sine 


-3 


or  for  a 3°  rotation  exi“^22“  -1. 37.10”  per  unit  time.  A consistent 


way  around  this  problem  is  to  calculate  the  coordinates  and  veloc- 
ities^of  th^  points  in  polar  coordinates  using  as  the  direction  of  the 


unit  r and  9 vectors  the  initial  polar  coordinates  of  the  points.  In 
particular,  if  the  initial  and  final  positions  of  a point  are  given 
by  the  polar  coordinates  (r,0)  and  (r',6')  with  0 measured  clockwise 
from  the  y axis,  then  the  x and  y velocities  are  given  by 


u^  = (r'-r)  sine  + r[e’-0]  cos0 


Uy,  = (r'-r)  COS0  - r[0'-0]  sin0 


Using  these  velocity  estimates  in  equation  1,  we  would  obtain  (using 
Figure  2)  for  a clockwise  rotation  of  0 per  unit  time 


"ll  "22  " "12  ° 


w = [0'-e] 


For  the  general  case  of  many  points  utilizing  least-squares  equations 
based  on  Equations  1,  it  can  be  shown  (Hibler  et  al.  , 197*^)  that 
adding  a constant  rotation  to  all  angles  only  changes  the  vorticity. 


Besides  the  special  example  considered  here,  this  rotation 
effect  can  be  illustrated  by  actual  data.  Such  an  Illustration  is 
given  in  Figure  3,  where  various  rotational  errors  in  the  coordinate 
transfer  of  five  common  ice  features  between  two  LANDSAT  images  was 
purposely  added.  As  can  be  seen,  in  agreement  with  the  simple 
example  of  Figure  2,  rotations  of  3°  can  cause  significant  undesir- 
able effects. 
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Figure  2.  EFFECT  OF  SPUBIOUS  ROTATIONS  ON  ESTIMATED  VELOCITIES 
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East-West  Strain  Rote 


North- South  Strain  Rote 
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Polor  Colculofion 
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Spurious  Rototion 


Figure 

3.  EFFECT  OF 

SPURIOUS  : 

ROTATIONS  ON  ESTIMATED 

STRAINS 

TABLE 

IV.  APPARENT 

STRAINS 

ON  LAND  MASSES. 

Units:  Percent  per 

day 

Photo  Pair 

e 

XX 

e 

X2. 

Shear 

Vorticity 

1 

-0.069  + 

0.106 

-0.276  + 

0.087 

O.Oltit  + 0.069 

0.063  1 0.069 

2 

-0.178  + 

o.oltl 

-0.033  + 

0.108 

0.2lt7  +.  0.058 

-0.386  + 0.058 

3 

0.02U  + 

0.0lt9 

0.098  + 

0.053 

-0.2lt6  + 0.036 

0.1lt3  + 0.036 

It 

-0.2U1  + 

o.iolt 

0.118  + 

0.098 

0.108  + 0.072 

-0.051  + 0.072 

5 

-0.11i2  + 

0.0U6 

0.159  + 

0.053 

0.109  * 0.035 

0.29*1  + 0.035 

6 

-0.020  + 

0.15T 

0.1t02  + 

0.089 

-0.078  + 0.090 

0.*25  + 0.090 

RMS 

0.139 

0.219 

0.159 

0.215 
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In  using  this  "polar"  technique  it  should  he  noted  that  as  long 
as  one  chooses  one  of  the  points  as  the  origin,  where  the  center  of 
rotation  is  located  has  no  effect  on  the  results.  This  is  because  all 
lines  connecting  points  will  rotate  by  the  same  amount  regardless  of 
where  the  center  of  rotation  is  located.  This  will  not,  however,  be 
true  if  one  chooses  for  the  origin  some  arbitrary  point  which  may  not 
actually  represent  a point  on  the  ice. 

Although  the  above  procedure  allows  strain  estimates  to  be  inde- 
pendent of  coordinate  transfer  errors,  there  will  nevertheless  still 
be  some  strain  errors  because  of  the  registration  error  in  maurking 
common  features  and  the  non-linearities  in  the  images.  Besides  being 
dependent  on  such  factors,  vorticity  estimates  additionally  depend  on 
errors  in  the  orientation  estimates  of  the  coordinate  systems  on  the  two 
images. 

To  make  a direct  estimate  of  these  errors  we  calculated  least- 
squares  strain  and  vorticity  for  the  photo  pairs  used  in  the  above 
transferral  study.  The  results  are  shown  in  Table  IV  where  we  list  the 
average  strain  rates  and  vorticities  (which  should  be  zero  since  there 
has  been  no  motion)  together  with  the  inhomogeneity  errors.  As  can  be 
seen,  on  the  average  the  strains  are  of  the  order  of  0.1  to  Q,2%  per  day 
which  would  represent  a 0.2  km  variation  over  100  km,  the  approximate 
size  of  the  arrays  we  were  using. 

However,  not  all  the  error  can  be  attributed  to  random  errors 
such  as  the  registration  errors.  If  this  were  the  case  then  the 
average  strain  rates  in  Table  TV  woxild  cluster  around  zero  with  the 
error  being  somewhat  larger  than  the  average  strain.  A more  consis- 
tent way  to  interpret  the  results  in  Table  TV  would  be  to  say  that  in 
addition  to  registration  errors  there  are  certain  distortions  in  the 
photos.  Assuming  such  distortions  are  approximately  linear,  the 
estimated  strain  would  be  equal  to  the  average  distortion  with  the 
residual  error  indicative  of  the  registration  error.  For  the  six 
photos  analyzed  the  average  residual  error  is  86  m which  is  almost 
exactly  the  stated  resolution  (80  m)  of  the  LANDSAT  multispectral 
scanning  system.  The  rest  of  the  strain  can  be  accounted  for  by 
saying  that  distortions  induce  about  100  m errors  per  100  km.  Conse- 
quently as  a practical  manner  for  estimating  experimental  errors  in 
strain  we  can  assign  equal  errors  of  about  100  m to  the  registration 
of  points  (which  is  in  agreement  with  Colvocoresses  (197**)  estimates 
of  geometric  fidelity)  and  to  image  nonlinearities  over  distances  of 
100  km. 

With  regard  to  vorticity  errors  the  rms  vorticity  in  Table  IV  was 
only  a little  larger  than  the  rms  strain  rate  indicating  that  rotation 
effects  were  relatively  minor.  However,  comparisons  of  the  angles  of 
the  line  of  longitude  passing  through  the  center  of  a lANDSAT  image,  as 
obtained  from  different  lines  longitude  marked  on  the  image,  indicate 
that  the  estimated  andles  often  differ  by  0.01  radian.  This  effect  is 
somewhat  ameliorated  by  taking  the  average  estimated  angle  from  two 
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longitude  lines  located  as  near  to  the  image  center  as  the  image  tic 
marks  allow.  Still,  hecause  of  such  errors,  it  would  seem  safer  to 
assume  the  experimental  vorticity  error  to  he  of  the  order  of  0.3%  per 
day. 


To  illustrate  the  application  of  these  techniques  to  actual  sea 
ice  drift  and  deformation  measurements,  we  examined  three  sequences  of 
images  from  March  1973  in  the  region  78°N' , 170°W.  The  starting  posi- 
tion and  day  for  each  sequence  is  shown  in  Figure  U.  Seq.  11  consis- 
ted of  7 imeiges;  Seq.  12,  6 images;  and  Seq.  13,  ^ images.  Four  to 
six  common  ice  features  were  digitized  on  each  sequential  image  pair 
and  the  average  drift  and  deformation  calculated.  Hie  drift  rates 
are  shown  in  Figure  5.  The  westward  drift  rates  generally  decrease 
as  one  moves  to  the  north,  a fact  which  is  in  general  eigreement  with 
the  clockwise  motion  of  the  Pacific  Gyre.  Also  the  drift  rates  are 
on  the  average  several  kilometers  a day  which  is  larger  than  the 
observed  transferral  error  of  about  1 km  per  day.  In  general,  these 
drift  rates  also  agree  with  average  drift  rates  of  ice  islsuids  which 
are  found  to  typically  drift  about  3.6  km/day  in  this  region  (Dunbar 
and  Wittmann,  1962), 


Figure  1*.  LOCATION  OF  INITIAL  IMAGES 
IN  SEA  ICE  IMAGE  SEQUENCES 
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Figure  5.  DRIFT  RATES  FROM  SEA  ICE  IMAGE  SEQUENCES 


DEFORMATION  RATE  ESTIMATES  FROM  AERIAL  IMAGERY 

Because  of  its  regular  orbit,  LANDSAT  produces  imsigery  which  has 
a consistent  orientation  and  linearity  that  makes  it  quite  useful  for 
estimating  deformation  rates.  Aerial  imagery  unless  corrected  by 
photogrammetric  techniques  will  have  orientation  differences  due  to 
changes  in  both  the  altitude  and  orientation  of  the  aircraft  during 
the  reconnaissance  run.  To  obtain  some  estimates  of  how  well  such 
luicorrected  imagery  can  be  used  to  meas\ire  distance  changes  bewteen 
ice  features,  and  hence  strain,  we  have  utilized  a series  of  aerial 
photos  take  at  91^**  m (30,000  ft)  by  the  NASA  Convair  990  in  April 
1975  over  the  AIDJEX  main  base  camp.  On  the  ice  two  i(0'  by  hO’  black 
cloth  targets  were  constructed  and  placed  approximately  7 km  apart . 
The  changes  in  the  distance  between  these  two  targets  were  measured 
to  an  accuracy  of  + 2 meters  using  two  theodolites  separated  on  a 
baseline  of  ss.1.5  km  in  length. 

Using  six  sequential  images  over  the  camp  relative  distance 
changes  between  the  targets  were  estimated  by  comparing  the  target 
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distances  to  distances  tietween  fixed  "immovable"  points  on  undeformed 
areas  of  a multiyear  ice  floe.  In  particular,  3 non  parallel  lines 
several  kilometers  in  length  were  identified  on  a given  floe.  The 
average  ratio  of  the  lengths  of  these  lines  as  determined  from  one 
photo  in  the  sequence  to  their  lengths  as  determined  in  an  earlier 
photo  was  determined  and  used  as  a correction  factor  to  be  applied 
to  the  distance  between  the  visual  targets.  The  comparison  of  actual 
and  photo  estimated  distances  are  given  in  Figure  6,  where  we  have 
asstamed  the  target  distance  on  the  first  photo  to  be  equal  to  the 
actual  distance.  As  can  be  seen,  the  photo  distances  show  a spread 
of  about  100  m or  1,3%  at  7 km.  Althoxigh  this  is  a significant  error 
it  is  still  small  relative  to  the  actual  distance  change  between  the 
targets  which  was  several  hundred  meters  3%) . Another  estimate 

of  this  type  of  error  was  made  by  examining  the  variation  from  photo 
to  photo  of  the  ratios  of  the  three  fixed  calibration  lines  among 
themselves.  If  there  were  no  errors,  these  ratios  would  not  change, 
whereas  in  practice  the  standard  deviation  over  6 photographs  varied 
from  1.7^  to  2.2%. 

Consequently,  these  limited  results  suggest  that  2%  is  a reason- 
able estimate  for  the  error  in  measuring  distance  from  high  quality 
aerial  imagery  without  photogrammetric  corrections.  Because  actual 
sea  ice  distance  changes  are  often  10^  when  measured  on  the  small 
scale,  such  error  can  occasionally  be  neglected.  An  example  of  such 
a case  would  be  a study  where  only  the  widths  of  the  leads  are  impor- 
tant. 
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CONCLUSIONS 


We  believe  that  this  study  has  shown  that,  when  analyzed  with 
appropriate  techniques,  LANDSAT  imagery  is  adequate  for  estimating  sea 
ice  drift  rates  at  locations  far  from  land.  The  transferral  geometry 
described  in  this  paper  gives  a simple  and  rapid  way  to  transfer  co- 
ordinates from  one  image  to  another  with  the  errors  in  the  transferral 
equations  being  negligible  over  the  transferral  distances  used.  In 
general,  the  errors  in  the  transferral  process  are  dominated  by  the 
deviation  of  the  actual  image  centers  from  their  stated  centers.  The 
effect  of  such  uncertainties,  however,  is  not  enough  to  mask  effects 
caused  by  typical  sea  ice  drift. 

As  regards  deformation  rates,  typical  results  show  that  it  is 
important  to  use  a least  squares  procedure  that  is  independent  of 
rotation  effects.  Also,  the  magnitude  of  errors  in  strain  rates,  due  to 
experimental  errors  and  non-linearities  in  the  ice  velocity  field, 
suggest  that  two  day  averages  of  strain  rates  should  be  used  when  esti- 
mating deformation  rates  from  LANDSAT  images.  When  estimating  deforma- 
tion rates  from  uncorrected  aerial  imagery  of  sea  ice,  typical  errors 
appear  to  be  of  the  order  of  2%  which  although  large  are  again  less  than 
observed  differential  distance  changes  as  measured  on  a small  scale 
( several  kilometers ) . 
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